Impact of a complex singlet: electroweak baryogenesis and dark matter 
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With the assistance of a complex singlet, and an effective operator involving CP violations, the 
dark matter relic abundance and baryon asymmetry of the universe have been addressed simul- 
taneously. We studied the electroweak baryogenesis mechanism systematically. The electroweak 
phase transition analysis indicates that the strong first order phase transition takes place by one- 
step or two-step type due to the dynamics of the energy gap between the electroweak vacuum and 
the vacuum of the complex singlet. The relation between the magnitude of baryon asymmetry of 
the universe and the phase transition type and strength has been explored in the framework of 


electroweak baryogenesis. 
I. INTRODUCTION 


The Standard Model (SM) has passed most experi- 
mental tests during the last 40 years. And the discovery 
of a 125 GeV SM-like Higgs particle by the ATLAS 
and CMS [2] collaborations at the LHC seems to provide 
the last missing piece of the SM. However, it has long 
been known that the SM has two obvious shortcomings, 
i.e., the explanations of the observed dark matter abun- 
dance and matter-anti-matter asymmetry of the universe. 


Firstly, the observed dark matter abundance is one 
issue that couldn’t be addressed in the SM. The exis- 
tence of dark matter (DM) has already been established 
by the observations of galaxy rotation curves and analy- 
sis of cosmic microwave background (CMB) etc. And the 
PLANCK [4] and WMAP [3] predict that DM constitutes 
about 26.5% of our Universe. Secondly, the matter-anti- 
matter asymmetry of the universe, i.e., baryon asymme- 
try in the Universe (BAU) as a baryon to entropy ra- 


tio 4 5] 
7 ~ (0.7 — 0.9) x 10729 , (1) 


couldn’t be addressed and predicted in the framework 
of the SM. The dynamical generation of BAU requires 
three necessary ingredients [6]: (1) violation of baryon 
number; (2) violation of both C- and CP-symmetry; and 
(3) departure from equilibrium dynamics or CPT vio- 
lation. Though the SM contains all these requirements, 
and the electroweak phase transition (EWPT) provides a 
natural mechanism for baryogenesis [6], the SM is unable 
to solve the problem for the reasons being listed below: 
Although the baryon number violation is provided in the 
SM by weak sphalerons [7}9], the departure from ther- 
mal equilibrium, provided by a strong first order phase 
transition, proceeding by bubble nucleation which makes 
electroweak symmetry breaking (EWSB), does not oc- 
cur in the SM [I0]. In the SM, no first order phase 


transition occurs for Higgs mass larger than about 80 
GeV [IHI3], which is far below the experimental bound 
of mp > 114 GeV from LEP [I4] and m),~125 GeV 
from the LHC [I] 2]. In addition, the CP violation in the 
CKM matrix is too small to produce a sufficiently large 
baryon number. New CP violations beyond the SM are 
required [I5]. 

This work is aimed at accommodating the observed 
DM relic density and the BAU simultaneously. The ob- 
served DM relic density reported by CMB anisotropy 
probes suggests weakly interacting massive particle 
(WIMP) [20)[21] and a feasible candidate for DM requires 
the extension of the SM. Based on CxSM [22], we con- 
sider one simple extension of the SM with a complex sin- 
glet (S) being supplemented. The singlet transforms triv- 
ially under the SM gauge groups and the imaginary part 
takes the responsibility of being DM. One very attrac- 
tive mechanism, i.e., electroweak baryogenesis (EWBG), 
wherein baryon number generation is driven by the CP 
asymmetry at the time of the EWPT!, is adopt. We 
take advantage of the complex singlet S to obtain the 
CP-violation required by EWBG through introducing a 
dimension-6 operator. Ref. [48] considered the analogous 
dimension-5 operator involving S/A, and Ref. [49] use 
S?/A? because of the Z symmetry S — —S needed to 
prevent decay of S, as befits a dark matter candidate. In 
our case, we consider the complex scalar S, thus we have 
SS1/A?, which ensures that the dark matter candidate 
preserves Zə symmetry and the CP violation source sur- 
vives after symmetry breaking. The relevant Lagrangian 
takes the form of 

QLH (1 +4 (a + ib) 


ie ss" tr +h.c. (2) 


1 Which should be the strong first-order EWPT to protect the 
generated baryon asymmetry from washout and suppress the 
sphaleron after the process is ended. 


where a, b are real parameters and A is a new physics 
scale. During the EWPT, the top quark mass gets a 
spatially-varying complex phase along the bubble wall 
profile, which provides the source of CP violation needed 
to generate the baryon asymmetry. Precision tests, es- 
pecially the electric dipole moment (EDM) searches, can 
probe directly CP violation relevant to EWBG. Using the 
polar molecule thorium monoxide (ThO), the ACME col- 
laboration reported an upper limit on the electron EDM 
(eEDM) recently [16], at 90% confidence level, an order 
of magnitude stronger than the previous best limit, 


Ide] < 8.7 x 10779ecm (3) 
This limit severely constrains the allowed magnitude of 
CP-phases in the Higgs couplings via Barr-Zee 
diagrams. We would like to mention that, one must take 
into account the tension between the CP-phase being re- 
quired to successfully implement EWBG and constraints 
from the eEDM. The preliminary results of the work are 
listed as follows: 


e The EWPT have been explored in two parameter 
spaces: triple couplings (c2, 61) and quartic cou- 
plings (dz, 62), as well as cubic and quartic cou- 
plings (2,c2). The behaviors of the strength of 
EWPT (ve/Te) and the energy gap AV ( the en- 
ergy difference between the Electroweak vacuum 
and the vacuum of the complex singlet) are found 
to be opposite: a smaller AV corresponds to a 
bigger ve/Te. Two types of phase transition, i.e., 
one- and two-step are studied, and the two-step 
phase transition always gives rise to a bigger ve/Te, 
which is found to the same as the real singlet case 
as explored in [23]. The dynamics of the phase 
transition could be characterized by AV to some 
extent. 


The BAU during the EWPT in the model has been 
explored. The behaviors of BAU as functions of 
triple couplings and quartic couplings match well 
with that of EWPT. The magnitude of BAU is 
proportional to the strength of the strong first or- 
der electroweak phase transition (SFOEWPT) for 
the one-step EWPT, and could be a little higher 
in the two-step EWPT case with the sign-flip be- 
havior appears at some larger ve/Te. 


The imaginary part of the complex singlet serves 
as the DM candidate. The smaller dark matter 
mass ( ma) region is severely constrained by the 
direct detection experiment LUX. A larger m4 can 
give rise to a relatively higher magnitude of the 
relic density. The EWPT in the higher magnitude 
of the dark matter mass region could be indepen- 
dent of ma. A benchmark scenario which accom- 
modates DM and EWBG is presented. In the sce- 
nario, the magnitude of the BAU could fall into 
the observed range of PLANCK [4] in the two-step 


SFOEWPT situation, and the CP violation phase 
needed for the EWBG is allowed by the constraints 
of ACME. 


Il. THE MODEL 


With the notation as in [22], the tree-level potential 
of the model is given by, 
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V (H,S) = sm? HH + 2 (HYH)? + 2 mi ISP 
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+(7 breidg? 4 2E gig)? +c.c.) , (4) 
where H and S are the SU(2) doublet and complex sin- 
glet fields. For simplicity, phases are chosen as: œ», = 0 
(T), ds, = 0 (r), and e = 0 (a). To get the min- 
imization conditions of the potential, it is convenient 
to represent the SU(2) doublet and complex singlet as 

= (0, h/V2) and S = (S +iA)/V2. Thus Eq. (4) 
recasts the form of, 


Aa, v2 
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Vo (h 8,4) = FR + oh 
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+I (bg + bı) A? + 4 1 tbs = bı) S? + pes 


x (9 + A?) + 924? + 7 (S* + A*) 


in which 6, and c2 can be both positive or negative. As 
could be seen from Eq. (5), the breaking of the global 
U(1) induces one Z2 symmetry for the imaginary part of 
S, i.e. A, thus makes a feasible DM candidate. Mean- 
while, the cubic and quartic interactions that are analo- 
gous to [23] are both kept for the EWPT study. 


The three minimization conditions of the potential, 


(OVo/Oh)|h=v,S=2,A=0 = 0, 
(OVo/OS)|h=v,S=2,4=0 = 0, 


(OVo/OA)|n=v,s=2,A=0 = 9, (6) 


with v = 246.2 GeV and « being the vacuum expectation 
values (VEVs) of the h and S, imply that, 


m? = vi2b1 02 2 aro 
2 2 2 (7) 
bə = by v261 vu /x v2e2 x d2 x ô2 v? 
4 2 2 2 
At the minima, the mass matrix of h and S is, 
M= Ew? taxu + v2 yj 
Eôru + v2 y6, idr? H Bae v2? f 
(8) 


o (S? ae A?) 


9 (5) 


and the DM candidate A, with mass m4 = bı — 
oor = VZ? /z, does not mix with h or S. The 
non-zero mŽg induces mixing between the SM Higgs 
and the real component of the singlet (S). The mass 
square of the “Higgs-like” particle mj, is given by: 
5 (mz +m2+/(m? — m2)? + 4(m? ,)?) for m? > mz, and 
a (mp, + ms — y (m — mz)? + A(m7,,)?) for mi < ms . 
For m? > m2, the “Higgs-like” eigenstate hı will be the 
heavier one, and the “Singlet-like” eigenstate hə will be 
heavier at the time of m? < m?. Demanding mp, = 125 
GeV, A could be given as a function of (c2, 61) or (d2, 42), 
as shown in Fig. |1|. The mass eigenstates h,,h2 are 
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FIG. 1: Top: A as a function of 6, and c2 from mp, = 125 GeV 
with relevant parameters being set as: x = 200 GeV, m4 = 
300 GeV, d2 = 1.4, 62 = 0.2; Bottom: A as a function of d2 
and d2 with mp, = 125GeV, relevant parameters are set as: 
x = 200 GeV, ma = 300 GeV, 61 = —350 GeV, and c2 = 100 
GeV. 


related with gauge eigenstates h and S as, 


eit. © 


The mass eigenstates h;,h2 couple to the fermions and 
gauge bosons via SM Higgs couplings reduced by a factor 
of cos@ and —sin@ respectively. We take the mixing 
angle to be —1/4 < ġ < 7/4 so that hı(h2) is always 
the “Higgs-like” (“singlet-like”) eigenstate. Using this 
convention, the mass matrix elements of Eq. could be 
expressed as, 


mi, = cos ¢*m;,, + sin om}, , (10) 
m? = sin om, + cos PM, ; (11) 
Mig = cosPsin d(mp,, — mi,,) , (12) 


which allows us to replace A, d2 and 6; with mp,,mn, 
and ¢. As for the vacuum stability condition, parameters 
living in d2 < 0 and 62 < —VAd2 will be excluded. 
For more studies on the stability problem in this kind of 
model, we refer to [2426]. 


Ill. EWBG ANALYSIS 


Before the EWBG is triggered, the early Universe is 
assumed to be hot and radiation-dominated, containing 
zero net baryon charge, and the SU(2); x U(1)y elec- 
troweak symmetry is manifest [33}36]. With the uni- 
verse cooling to temperature below the electroweak scale, 
the Higgs field acquires a non-zero expectation value and 
thus spontaneously breaks the SU (2); xU (1)y symmetry 
down to U(1)zm™. During this phase transition process, 
the EWBG takes place. 


Bubble (of the broken phase) formation and growth 
begin at nucleation temperature Ty( < To, the critical 
temperature of SFOEWPT). Baryon creation in EWBG 
takes place in the front (and vicinity) of the expanding 
bubble walls [37]: With CP violation implemented, par- 
ticles in the plasma (symmetric phase) scatter with the 
bubble walls and generate CP (and C) asymmetries in 
particle number densities; and these asymmetries diffuse 
into the symmetric phase ahead of the bubble wall and 
they bias electroweak sphaleron transitions [I5] to pro- 
duce more baryons than antibaryons; and some of the 
net baryon charge being created is then captured by the 
expanding wall into the broken phase; at last, inside the 
bubbles(broken phase), the rate of sphaleron transitions, 
which provides (B + L)-violating processes and could 
wash out the baryons created, should be strongly sup- 
pressed, i.e., the SFOEWPT is needed. 


A. SFOEWPT realization 


The dynamic of the EWPT is governed by the effec- 
tive potential at finite temperature, i.e., Vesp (h, S, A, T) 
in our model. After minimizing Vess (h, S, A, T), the 
VEVs of three scalar fields h, S, and A at each tem- 
perature could be obtained. When a first order phase 
transition takes place, a local minimum with (h) 4 0 
develops and becomes degenerate with the electroweak 
symmetric minimum as the temperature decreases to the 
critical temperature Te, and the two minima are sepa- 
rated by an energy barrier. To protect the BAU gener- 
ated by EWBG from being washed out, the strong first 
order condition [17], 


> 1. (13) 


Ve 
Te 


with ve being the Higgs VEV of the symmetry broken 
minimum, is necessary. Here, when one uses the condi- 
tion, the following three most important theoretical un- 
certainties need to be kept in mind [18]: ve in the condi- 
tion is indeed the value of the classic field at minima at 
T.(¢(T.)) which depends on the gauge choice; the gauge 
dependence of Te does not compensate that of ¢(Te); and 
the choice of the washout factor is subject to additional 
uncertainties. 


The standard analysis of Vess (h, S, A,T) involves 
tree level scalar potential (Vo (h, S, A)) and Coleman- 
Weinberg one-loop effective potential (Vi (h, S, A)) at 
T = 0, the counter terms V., being chosen to main- 
tain the tree level relations of the parameters in the 
Vo (h, S, A), and the leading thermal corrections being 
denoted by V (h, S, A, T), 


Veff (h, S, A, T) = Vo (h, S, A) +V (h, S, A) +V (h, S, A,T) 


Here, Vo is the same as the tree level potential Eq. (5), 
Vı is the Coleman-Weinberg one-loop effective potential 
at T = 0 and could be expressed in terms of the field- 
dependent masses m; (h, s, A): 


mi (h, S, A m? (h, S, A 
V (h, 5, A) = err ) im ( D 2) - c| 
(15) 
in MS scheme and Landau gauge. Here, the sum i 
runs over the scalar, fermion, and vector boson contri- 
butions. n; is the number of degrees of freedom for the 
i-th particle, with a minus sign for fermionic particles. 
Q is a renormalization scale which we fix to Q = v. 
We take C; = 4 (C; = 2) for the transverse (longitu- 
dinal) polarizations of gauge bosons, and C; = 3 for all 
other particles. The field-dependent masses are given in 
the Appendix. B. With V; being entailed in the poten- 
tial, the minimization conditions Eq. (6) will be shifted 
slightly, and the relations Eq. (7) do not hold as well. To 


maintain these relations, the counter terms Vet should be 
added [38} [39], and cast the form of, 
Vee = Om@h? + MS? + MZA? , (16) 


in which, the relevant coefficients are determined by, 


1 OV, 

bmi = — 55 Gp InavSae,A=0 ; (17) 
1 ov, 

ôm = -— > =v,S=2,A=0 > 18 

Mog on OS [h= „S=x,A=0 (18) 
1 8? Vi 

m3 = -3 ggz lh=v,S=2,A=0 ` (19) 


Thus, the VEVs of h and S as well as the dark matter 
mass will not be shifted. We do not include more com- 
plicate terms to compensate the shift of mass matrix of 
S and h, because these shift effects are basically small. 


The finite temperature component V (h, S, A, T) is 
obtained as 


M2(h, A, S,T 
WASD) oop 


T4 
V (h, S,A,T) = a2 So niJB,F ( 


where the functions Jg p(y?) are given by 


Js rly) = f dx x° In [1 exp (- x? +y)| , (21) 


with the upper (lower) sign corresponds to bosonic 
(fermionic) contributions. And M?(h,S,A,T) is given 
in terms of the field dependent mass at T = 0 and the 
finite temperature mass function II; as, 


with II, being listed in Appendix. B. 


The behaviour of EWPT strongly depends on the 
shape of the potential at zero temperature. Fig. [2| shows 
contours of Vo in the h — S$ plane in two typical situa- 
tions. Since the Vo does not develop a minimum at A Æ 0 
when the dark matter mass is large as will be explored in 
Sec. [IV B| we conduct the analysis in the A = 0 plane. In 
the top panel, the electroweak vacuum is the only local 
minimum of the potential Vo. For this situation, the uni- 
verse might transform from the high-T symmetric phase 
to the electroweak symmetry broken phase directly. In 
the bottom panel, an additional local minimum besides 
the electroweak vacuum appears in the direction of S at 
h = 0. For this situation, there is a possibility that the 
universe transforms to this phase at the first step and 
then to the electroweak symmetry broken phase at the 
second step. And in fact, the first investigation of the 
real singlet case shows that the two-step transition 
(with the second step from the high temperature singlet 
vacum to the electroweak vacuum) gives rise to a higher 
magnitude of EWPT strength in comparison with the 
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FIG. 2: The contours of Vo (in units of GeV*) in h — S ( 
in units of GeV ) plane. With parameters being fixed as: 
xz = 200 GeV, 6:=-300 GeV, d2=1.4, d2=0.2, A=1.0, and 
c2 = —10 (100) GeV for the top (bottom) panel. 


one-step transition case. ? Since the potential at this 


additional minimum (in the direction of S') can be very 
close to that of the electroweak vacuum, the second step 
can happen at a very low temperature, therefore a strong 
first order phase transition could be realized easily. Of 
course, the phase transition could also be the one-step 
case, as long as the energy gap between the minima in 
the direction of S and h is too large to be overcome by 
the T dependent thermal correction at high temperature. 
For the similar analysis in a triplet scenario and Next-to- 
Minimal Supersymmetric Standard Model (NMSSM) we 


2 It is also feasible to study fermion dark matter together with 
strong first order phase transition {19}. 


refer to [41] and [42}{44). 
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FIG. 3: EWPT with parameters being fixed as: x = 
200 GeV, ma = 300 GeV, 62 = 0.2, d2 = 1.4, and À be- 
ing solved from ma, = 125 GeV as shown in Fig. [1] The blue 
shade is the regions where mpi = 125 GeV couldn’t be ob- 
tained. The grey region depicts that the electroweak vacuum 
is a local but not global minimum. The light orange and red 
regions are parameters space where one-step and two-step first 
order phase transitions take place. The second order phase 
transition happens in the light blue region. Top: The green 
line is to divide the parameter spaces into the regions where 
the minimum in the direction of S exists (above) and disap- 
pears (below). The black contours in the region above the 
green line are the contours of AV = Vs — Vh in units of 10° 
GeV?; Bottom: The contours of ve/Te, cos@, and Mh are 
shown with orange(and red), pink and brown lines. The re- 
gion excluded by cos ¢ < 0.8, which characterizes the mixing 
between h and S, is covered by light grey shade. 


We take advantage of the CosmoTransitions pack- 
age [45] to investigate the phase transition numerically. 
Fig. [3] shows the result in the cp — 6, plane. In the top 
panel, the contour lines of Vo are plotted to help under- 
stand this pattern qualitatively: for parameter regions 
bellow the green line, the electroweak vacuum is the only 
minimum of Vo, thus only one-step phase transition could 
happen. The contour AV = 0.6x 108 GeV‘ separates the 
one-step from the two-step type phase transition. When 
the difference between the two minima is bigger than 
that value, the T dependent thermal corrections could 


only changes the magnitude of AV a bit thus only the 
one-step phase transition could be obtained. When the 
magnitude of AV gets smaller, the phase transition type 
changes from one-step to two-step. The bottom panel 
indicates that, with the increase of c2 and 61, one gets 
smaller v./T. for the one-step phase transition, and a 
second order phase transition is obtained at last. For the 
two-step cases, ve/Te gets larger with the increase of c2 
and 6; and the SFOEWPT condition ve/Te > 1 could 
always be satisfied. The pattern of the contours of ve/Te 
is opposite to that of AV, i.e., the smaller AV leads to 
the larger ve/Te. 


62 


FIG. 4: With parameters being fixed as x = 200 GeV, m4 = 
300 GeV, 61 = —350 GeV, and c2 = 100 GeV. The descrip- 
tions of the two pictures are the same as in Fig. 


Fig. |4| illustrates the dependence of the behavior of 
phase transition on the quartic couplings ô and dz. For 
the one-step phase transition, the magnitude of ve/Te 
becomes higher and higher with the decrease(increase) 
of 62 (dg). For the two-step type phase transition, the 
value of ve/Te increases with the two couplings getting 
larger. 


As depicted by the bottom panels of Fig. B]and Fig. [4] 
for our parameter choice, one-step type phase transition 
parameter regions are severely constrained by the mix- 
ing angle by fitting current Higgs measurements at the 


LHC [31I], while most two-step ones survive. And to 
make a decisive conclusion on this topic, more data are 
expected from the run II of the LHC as well as other 
future colliders, such as the International Linear Col- 
lider (ILC), TLEP (renamed as FCC-ee recently), China 
Electron Positron Collider (CEPC), or the Super proton- 
proton Collider (SPPC), etc. Here, we would like to men- 
tion that the electroweak precision observables are sup- 
posed to constrain the Higgs and singlet mixing angle ¢ 
and the heavy Higgs mass mp, [67], which would 
also bound the SFOEWPT and BAU favored regions. We 
plot the contours of the mixing angle and masses of the 
“Singlet-like” particle in our parameter regions to take 
into account this thing. For detailed analysis we refer 


to BI]. 


B. BAU gereratation 


To get a more comprehensive understanding of 
EWBG, we solve the transport equations as explored 
in Ref. [40] directly, in which the speed and profile of 
the walls are the core ingredients. With the nucleating 
bubbles reaching a sizable extent and expanding with a 
constant velocity, it is appropriate to boost into the rest 
frame of the bubble wall and proceed all calculations with 
a planar wall. With z being the coordinate transverse to 
the wall, h(z), S(z) and A(z) are the expectation values 
of h, S and A inside and outside the bubbles. To deter- 
mine them precisely, one needs to minimize the Euclidean 


action [46]: 
_ [T alton Lao 
J= J [Gen +5 (0.8) 
1 
+5(0:A)* + Vess (h, S, A, Te) : (23) 


As a good approximation, we estimate the bubble wall 
profiles in the form of, 


h(z) = sell —tanh(z/Lw)), (24) 
S(z) = at lAr + tanh(z/Lw)), (25) 
A(z) = ac + 5 Aac( +tanh(z/Lw)) . (26) 


Here, Ly is the width of the wall, ve, £e, and ae are 
the VEVs of the electroweak symmetry broken phase at 
T., and Ax, (Aac) is the total change of (S) ((A)). For 
the scenario of this work, A hardly gets VEV during the 
whole phase transition process, therefore, A(z) is zero, as 
will be explored more in section [IV B| The wall width is 
approximated with Ly ~ (Ax? + Aa? + v2) /8V,, ÆJ, 
in which Vz is the height of the potential barrier at Te. 
The bubble wall profiles are shown in Fig. 


As mentioned in Sec. [I| we introduced the CP violat- 
ing top quark mass term, and for definiteness we fix a = 0 


new phase old phase 


<h> 


FIG. 5: The wall profiles at Te with new (old) phase corre- 
sponding to Electroweak broken (unbroken) phase. 


to maximize the CP violation. The top quark acquires 
a spatially varying complex mass inside the bubble walls 
during the phase transition, being obtained as 


m(z) = BO (1 Po ai 2) = [miO . 
(27) 


following the method used in. 49]. We will verify if 
the nontrivial phase 6(z) is sufficient to source the baryon 
asymmetry in the following studies. The energy scale will 
be chosen as A = 1000 GeV in relevant calculations of 
this work, i.e., eEDM and BAU.? 


With the top quark profile of our model being ob- 
tained, we solve the transport equations ( see Apendix 
D for more details) to get the chemical potentials of 
left-handed SU(2) doublet top quark( puz,2), left-handed 
SU(2) doublet bottoms (u»,2), left-handed SU(2) singlet 
top quark (j1ic,2), and Higgs bosons (un,2), and the cor- 
responding plasma velocities. The chemical potential of 
left-handed quarks cast the form of, 


1 
a + 4K 15) Mo,2 — 2K1 tMee,2 , 
(28) 


1 
LB, = 51 +441 +) pe + 


where K factors are thermal averages [40]. 


Finally, The baryon asymmetry is obtained using 
405. ws _ 
“ep z dz MB, (z á , (29) 


ae An, 9T Vw Gal 
where Pws is the weak sphaleron rate and v = 
45. ws/(4Uw). The effective number of degrees of free- 
dom in the plasma is g* = 106.75. And we take the wall 
velocity to be vw = 0.1 in our calculation. 


1B = 


3 Here, we would like to mention that the energy scale A charac- 
terizes the size of 6(z) to some extent. And a larger A always 
lead to a smaller (z) and results in a lower magnitude of ng 
and de/e. The dimension six operator in Eq. induces negligi- 
ble effects in the EWPT process, and the wall profiles shown in 
Fig. [5]ensures that the operator’s contribution to m;:(z) is small 
in comparison with the renormalizable one 3 
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FIG. 6: The correlation between baryon asymmetry and the 
EWPT parameters. Parameters are fixed as ve = 150 GeV, 
an = 0, and Tc = 100 GeV(ze = 150 GeV, za = 0, and 
an = 100 GeV) for the top (middle) panel. The magnitude 
of ne as a function of x, and rp, ve and Te, and ve/Te are 
shown in top, middle, and bottom panels . 


To get a first glimpse of the behavior of the baryon 
asymmetry with respect to the parameters of EWPT ( 
i.e., the critical temperature T, and the field expectation 
values Uc,@¢,tp.--), we perform the calculation of BAU 
with these parameters being input directly. We need to 
note that the values of ve and Te used here are not phys- 
ical parameters we get during EWPT as explored in the 
previous section. The top panel in Fig. [6Jillustrates how 
the expectation values of S affect the magnitude of BAU, 
with x. and zp denoting the VEV of S at electroweak 
symmetry broken phase and symmetric phase. It depicts 
that ng gets larger with x, getting smaller for x. > 0, 
and the situation changes to be opposite when ze < 0. 
ng changes its sign near £a = £e. In compare with that 
of zn, Ng relies on x, more, especially when |x,| is small. 
The middle panel indicates how 7g depends on ve and 
Te. The dependence of 7g on ve is simple and obvious: a 
larger value of 7p corresponds to a larger ve as expected, 
since the source of CP violation depends partly upon the 
amount of variation of |m:|(~ ve) inside the bubble wall. 
The dependence of 7g on Te is more complicated. ng 
have a almost linear correlation with T, for larger Te, 
and the sign of 7g flips for smaller Te. The bottom panel 
gives the magnitude of BAU as a function of v./T., which 
indicates that some value of ve might induce unphysical 
results, i.e., ng < 0. Thus, a larger ve/Te may not leads 
to a bigger nz. 


Then, we use the results obtained in the section. 
as the input parameters and calculate the magnitude of 
the BAU being generated during SFOEWPT. The results 
are shown in Fig. |7| which depicts that for the one-step 
case the behaviors of 7g are very regular and indicates 
one good correspondence with the behavior of ve/Te, i.e., 
a simple linear correlation. While, a more complicated 
correlation for the two-step case exists. The sign of ng 
flips and flips again with the increase of ve/Te in the 
two-step case, which depicts a similar behavior as that of 
Fig. [6] 

As could be observed in the two panels of Fig. |7| a 
larger baryon asymmetry can be generated during the 
two-step type phase transition process. And, 7x 1071!! < 
ng < 9 x 1071! turns out to live in the two-step type 
EWPT favored region for the parameters set. Last but 
no least, we would like to go more deeper to the inner 
mechanism of the sign flip of np and relevant behaviors 
with respect to ve/Te as shown in Fig.|7} The BAU mag- 
nitude is highly related with that of up, (see Eq. 29), 
which is supposed to be dominated by the behavior of 
the source term S; + . To illustrate the physical picture 
more clear, we take some typical benchmark points be- 
ing used in Fig. [7| and plot the magnitude of upg, and 
the dominated source term S; with respect to the spa- 


4 Which is the same as that of the Sọ in [40], here we denote it as 
S+ to depict that it is the top source term. 
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FIG. 7: The BAU magnitude nz as a function of ve/Te, with 
orange(red) points depicting the one-step(two-step) phase 
transition cases. Parameters are set as: x = 200 GeV, ma = 
300 GeV, 62 = 0.2 ( 61 = —350 GeV) and d2 = 1.4 
(c2 = 100 GeV) for top (bottom) panel. . 


cial coordinate z, see Fig.|8| The oscillation amplitude of 
UB, and S; are found to be increasing with the increase 
of ve/Te. With the increase of the oscillation amplitude 
of ug, and St, ng increases at first and decreases lat- 
ter. And the volatility of the oscillation of ug, and Si 
give rise to unphysical BAU magnitude, i.e., the negative 
value of np appears (see the red line in Fig. [8). 


C. EDM 


In this section, we use the eEDM search to constrain 
the CP violation phase. The eEDM contribution is dom- 
inated by the Barr-Zee diagram [50], as illustrated in 
Fig. [9] Where, the two scalar mass eigenstates(h1,2) give 
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FIG. 8: Top and middle panels are the plots of we, and S; 
as a function of the coordinate transversing to the wall (z), 
with benchmark points being given in the bottom panel. 


two opposite sign contributions, obtained using ° 


de 16 a me 


e 3 (4r)3 v? 


oe Z [a (m/mi) -G (m?/mi) | , (30) 


x 


5 Here, only the dominate top quark contribution and hyy inter- 
action have been considered. As for the most general case, and 
possible cancellation mechanism we refer to ; 


h 1 / hy / 
/ 7 
> Z > > 


E e e 


FIG. 9: Barr-Zee Feynman diagram. 


where Me are the electron (top quark) masses, and Z 
characterizes the mixing between h and S [52], being ex- 
pressed as: 


Z = cos ġsin ġ = 


and the function G(z), 


ewe: f a, ab jo) (32) 


The magnitude of eEDM, in the parameter regions of 
Fig.|7| decreases with the increase of 61(d2) and is below 
the limit given by the experimental results of ACME [16], 
which means that the CP violation phase in the param- 
eter spaces is allowed by the eEDM experiment ACME. 


IV. DARK MATTER AND ELECTROWEAK 
BARYOGENESIS 


In this section, we investigate the possibility to ex- 
plain the dark matter relic abundance and EWBG in the 
same parameter spaces in the model. To achieve this 
purpose, we go through the dark matter issues at first. 
Then we explore the relation between the dark matter 
mass and the EWPT. The explanation of the the dark 
matter relic abundance and the EWBG will be acquired 
simultaneously at last. 


A. Dark matter 


In our model, the thermal relic density Q 4h? is mainly 
controlled by the pair annihilation cross section of A, i.e., 
Qah? ~ 1/(ov), in which (cv) is the thermal average 
of the product of annihilation cross section and relative 
velocity, with relevant Feynman diagrams being shown 
in Fig. 

Direct detection experiments are searching for dark 
matter scattering off atomic nuclei, and in the model we 
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FIG. 10: Annihilation channels of dark matter particle pairs. 


only have the Higgs mediated spin-independent scatter- 
ing cross section which constrains our parameter spaces. 


; : : 400 600 800 1000 1200 1400 
The scattering cross section of A with a proton for our 
model casts the form of, ma (GeV) 
ah? 
2 
ms, JAAh, CO8@  JAAh Sind 200 am 
Odd = z% 2 D | 
2rv? (mp =F ma) Mhi Ma, | 
2 | 
x (Jn + heat foo + = (fe)) | (33) 100 SA 012 | 
oe > | 0. A 
which is the same with Eq. (25) in [24] except that the 2 | 
couplings gAAh, and gaan, are different as shown in the = oF | 0. - 
Appendix. B. > | k 06 
Fig. [11] depicts the relic density as a function of m4 | 
and 62 (cg). The magnitude of the relic density in- -100% | | 
creases with the increase of m4, since the annihilation | 0.02 
cross section is inversely proportional to the dark matter | \ 
mass. The LUX experimental results mostly constrain —200 = i ; raci 
the smaller dark matter masses region, because the oqa 400 600 800 1000 1200 1400 
has the dark matter mass suppress effect, as could be ma (GeV) 
seen in Eq. The ô> and cz affect the magnitude of 


the relic density and direct detection of the dark matter 
mainly through the couplings gaan, and g9AAhs 


B. The dark matter mass and the EWPT 


In this section, we explore the relation between the 
dark matter mass and the EWPT. 


e The dark matter mass, T = 0 and T 4 0 vacuum 
structure 


The full analysis of the vacuum structure in our 
model is not easy, since there are many extreme 
points in the three dimensional space and the lack 
of Zə symmetry of S makes these extreme solutions 
very complicated. In the numerical scan of the 
parameter space in the EWPT part, we find that 
in most cases there are no local minimum at A Æ 0, 
whether at T = 0 or T Æ 0, especially when m4 
is large. This can be understood partly through 


FIG. 11: Dark matter relic density Qh? as a function of c2 
(62) and ma with x = 200 GeV, ¢ = 0.6, and Mr, = 250 GeV. 
The contours of 04h? are shown by purple lines with cz = 30 
GeV (62 = 0.2) for the top (bottom) panel. The parameter 
spaces excluded by LUX are shown by the orange region. 


the following arguments. Suppose there is a local 
minimum (vp,Us,U4), it should also be a minimum 
in the S = v,,h = vp subspace. In this subspace, 
the potential Vo(vp,,vs,A) has a local minimum 
when the coefficient of A? term is negative: 


b2 + bı 


V2c2 də 2 62 9 
! st } 4 
5) 6 UV 4 Us A v2 < 0. (3 ) 
Using the Eq. |7| we rewire the above condition as 
2 d ô 
mi V2 (x—vs)4 rac a?)+ z (eh v?) <0. (35) 


As a nature expectation, Vs ~ x and vp ~ v. So 
when m4 is large, the condition of Eq. [35] can be 
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easily violated and a local minimum at A 4 0 does 
not exist in this subspace, neither in the whole 
space. At T Æ 0, thermal contributions stabilize 
A field at the origin so A # 0 minimum does not 
exist at high temperatures. 


EWPT and the dark matter mass 
With A = 0, the Eq. [5] recasts the form of 


2 
2 
Vo(h, S) =4 h? + A ag 2 5185 21282 
1 V2 d2 
= (bə — 2 aa 3 M4 g4 
+7 (ba by) S* + eS t , (36) 


Thus, the vacua gap AV = Vp — Vs at tree level 
could be obtained as, 


my 1 ca(v3 — z?) — do(vt — zł) 
A = b 2 s s 
K A bles 6/2 16 
S 2 2,2 
bə bı re v 0041 Yr 69 AÀ v4 (37) 
4 4/2 8 16 


with vs being the solution of 0gVo(0, S) = 0. Since 
the dark matter mass could only affect the com- 
bination bı + b2, and the combination b2 — bı is 
a constant with other parameters being fixed (see 
Eq. F), one finds that the dark matter mass is ir- 
relevant to the vacuum energy gap, thus one could 
expect that the strength of EWPT (inversely pro- 
portional to AV) won’t be affected by which. 


Based on the above two arguments, we found that the 
EWPT could be independent of the dark matter mass in 
the relatively larger dark matter mass regions. 


C. EWBG and DM 


The higher magnitude of the relic density requires 
a relatively larger dark matter mass, about ma ~ 1000 
GeV. And in this mass region, as explored in the last sec- 
tion and the section [ITB] the EWPT and the generation 
of BAU is independent of m,. In order to accommodate 
the dark matter relic abundance and the EWBG simul- 
taneously we assume the dark matter mass ma = 1500 
GeV. 


We present in Fig. [12] the combined result of EWBG 
and DM relic density. All parameter regions in this 
plot are under the direct detection constraint from LUX. 
Fig. [12] indicates that the behaviors of one- and two-step 
phase transitions are dominated by 62 and c2 respec- 
tively. The dark matter relic density could be affected 
by 62 and c2, since they enter into the couplings of the 
dark matter particle and Higgs, i.e., ga Ahı (h2), Which de- 
termines the dark matter particle pairs annihilation pro- 
cesses. A larger value of c3 leads to a larger Q.4h?, and 
the correct dark matter relic density observed by [3] 4] 
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could be obtained in some parameter spaces of d2 and 
co. As for the strength of EWPT and the magnitude 
of BAU, we found that the two-step SFOEWPT case is 
more favored. The magnitude of eEDM is obtained to 
be |d-| = 4 x 10~?° ecm, below the upper limit given by 
ACME. 


C2 (GeV) 


FIG. 12: The contours of the relic density of the dark matter 
(Q4h7) and the magnitude of ve/Te with x = 200 GeV, ¢ = 
0.6, ma = 1500 GeV, and ma, = 250 GeV. The contours of 
Qah? are shown by purple lines and the light orange(red) lines 
depict the value of ve/Te. The green points are parameter 
spaces with appropriate BAU magnitude 0.7 x 107!° < ng < 
0.9 x 107". 


V. COMMENTS ON THE CALCULATION 
METHODS INVOLVING SFOEWPT AND BAU 


As pointed out by [18], the SFOEWPT criteria ve/Te 
calculated in the traditional mathod is subject to gauge 
problems, as the ve and Te both involve gauge depen- 
dence problem. The EWPT is performed in Landau 
gauge in this paper and one needs to take this into ac- 
count when interpreting our results. We have checked the 
robustness of the analysis in the model by following the 
gauge independent treatment introduced in for part 
of the parameter spaces favored by the EWBG and DM 
explanations in Fig. We find that the magnitude of 
Te drops a lot and the corresponding v(Te)/Te is higher 
than ve/Te for the two-step case. 


The numerical calculation of the BAU is based on the 
source terms and chemical potentials transport equations 
being computed within the WKB approximation. The 
results being given might be different to some extent in 
comparison with the results calculated using the quan- 
tum transport equations from the close time path (CTP) 
formulation, wherein the CP violation source is derived 
using the “(vacuum expectation value)vev-insertion” ap- 


proach [73H76]. 


VI. DISCUSSIONS AND CONCLUSION 


By supplementing a complex singlet to the SM, DM 
relic abundance and EWBG could be accommodated si- 
multaneously. Both two ingredients of EWBG, i.e., the 
SFOEWPT and new CP violation phase, have been ex- 
plored at the same time. The imaginary part of the com- 
plex singlet (A) makes the DM candidate due to the re- 
duced Z symmetry from a global U(1) symmetry. 


The phase transition could be one- and two- step 
types as found by the study of the energy gap AV. In or- 
der to reveal the effects of cubic and quartic terms on the 
EWPT, we performed EWPT analysis numerically using 
CosmoTransitions package in the parameter spaces of 
(c2, 61), (d2, 62) and (d2, c2). Both one- and two- steps 
first order phase transition are obtained, and the two- 
step case favors SFOEWPT more than that of one-step 
case. And this study is consistent with the analysis of 
AV, i.e., a lower magnitude of AV results in a larger 
strength of EWPT. 


The newly introduced CP violation makes the EWBG 
mechanism viable, and the CP violation phase is found 
to be allowed by constraints of ACME. The general and 
specific calculations of ng shows that the magnitude of 
BAU is highly related with v and T,, and the increase of 
ve/Te does not always cause the increase of ng, because 
ve/Te affects the CP violation source. 


The magnitude of the relic density of the DM is found 
to increase with the increase of the dark matter mass. 
The smaller DM mass regions are severely constrained by 
direct detection experimental results. And the EWPT is 
irrelevant to the DM mass in the larger DM mass regions. 
At last, we present a scenario in which both correct DM 
relic density and successful EWBG could be achieved 
in some parameter spaces of d2 and co considering the 
eEDM experiment ACME. 


There is a possibility that baryogenesis and dark mat- 
ter are generated at the time of EWPT simultaneously 
with the dark matter sector being asymmetric dark mat- 
ter , thus shed light on the path to explain the ra- 
tio between the observed abundance of DM and baryons 
in the present day: Qpy/Qpg ~ 5.4 [4]. The key in- 
gredient for the test of this kind of model is the colliders 
search of the Higgs triple coupling and the 
introduced CP violation [62H66], that will be postponed 
to the future works. 


Appendix A: Field dependent masses 
Field dependent masses matrix is given by, 


mi s a = eigenvalues (M°) , (A1) 
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with 
m? (h, S, A) mẹ?g(h,5S,A) m2, (h, S, A) 
M? = | m?s (h, S, A) Mg (h, S, A) m34 (h, S, A) )| , 
mia (h, 5, A) m34 (h, 5, A) m4, (h, S, A) 
and matrix entries are given by 
m? V2ô1 ô2 3A 
May = z tg Stg (S? +A?) + Oh? (A2) 
mg = =) + 2s, (A3) 
bi +b2 d V2c2 3d2 
2 os, 2 
Mog = 5 av A* 4 6 S4 1 S 
e (A4) 
2 b, 
2 d 
ma, = Tes, + SAS, (A6) 
ee nth +e | "a ny 
42242 , (A7) 
4 
and other field dependent masses given by, 
m? V2ô1 ô2 À 
m2, = ar S4 7 (S? 4 A?) + Zh? (A8) 
1 
m? = =yh? ’ (A9) 
1 
mys = 39h, (A10) 
1 
mz = 5 (9° +9") h (A11) 


with G being goldstone fields. 


Appendix B: Scalar triple and quartic couplings 


—2i(őz2v cos 6/4 + (V2c2/12 + dox/4) sing) , 


JAAh, = 
Ihihihy = —6ilAv cos? 6/4 + (dox /4 + V2c2/12)) sin? o 
+(d9@/4 + V261 /8) cos? dsin ¢ 
+ô2v cos osin? 6/4) , 
JAAhy = —2i(—dgusin ¢/4 + (V2c2/12 + d2£/4) cos 6) , 
Jhihiha = —2i(—3Av cos? dsin 6/4 + 3sin ¢? cos o(d2£/4 
+V2c2/12) — (5ox/2 + V251/4) cos ¢ sin? b 
+62u/2sin ¢ cos? $) , 
JAAhth, = —4i(d2 cos” /8 + dz sin? ¢/8) , 
Jhahoha = —6i(—Avsin? 6/4 + (dox/4 + V2c2/12)) cos? ġ 
+(692/4 + V281 /8) cos p sin? ¢ 
—dgu sin ġ cos? ¢/4) , 
JAAhshy = —4i(d2 sin? 6/8 + dz cos” 6/8) , 
Ghohohy = —2i(3 cos? dsin o(doa/4 + V2c2/12) 
+3 sin? ¢ cos @du/4 — (2x /2 + V26,/4) 
x cos” dsin ¢ — ôzv sin? ġ cos $/2) , 
JAAhihə = —2i(sin dcos d(d2g — d2)/4) . (B1) 


Appendix C: Thermal masses used in EWPT 
analysis 


Thermal masses of bosonic field relevant to EWPT 
study are listed bellow, 


Ig = (397/16 + g’?/16 + A/8 + y?/4 + 52/24)T? , 
Tan = (397/16 + 9/16 + d/8 + y? /4 + 62/24)T? , 
Iss = (do/12 + 62/48)T? , 

aa = (d2/12 + 6/48)T? , 

Mk = I} = 119/67" 


m7, =m =0. (C1) 


Appendix D: Transport equations 


For the case when the top quark acquires a com- 
plex mass during the electroweak phase transition, m = 
Im(z)l]et? C), the source terms and transport equations 
were derived in Ref. [40] using the WKB approximation . 
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Here we briefly summarize the approach with notations 
adopted in Ref. [40]. With vu denoting the wall velocity, 
the plasma velocity u and the coefficients K;, K; are de- 
fined by thermal averages. And the chemical potentials 
are expanded to second order in gradients. 

The source terms induced by the top quark are 


Sy = K70'm7 py, 
So = —VyKg(m?6') + vy Kg m? (mY , 
Su = —Kipm?6'ul, . (D1) 


in which the primes denotes derivatives with respect to z 
(the direction perpendicular to the wall) and the CP-even 
first order perturbations u1, 4, can be solved by 


Vy Ka, + Up Ko(m?)' py + uy = Tt" ua = Uw K3 (m?) ’ 
—Kypi + Vw K5ul, + VwK6(m?)'uy +T%ty, =0. 


And, the Sọ is the dominated one, effects induced by Su 
and S,, are negligible. 


With the weak sphaleron rate [69], the strong 
sphaleron rate [70], the top Yukawa rate [68], the top 
helicity flip rate, the Higgs number violating rate [68], 
the W scattering rate, the quark diffusion constant 
and the Higgs diffusion constant being given by, 


Tus = 1.0 x 10-®T, T,, = 4.9 x 10747, 


m?2(z, T) 

T, = 4.2 x 10-7 Tm = 

y a 63T ’ 

2 
= miy (z,T) __ tot 
Tp, = 50T o] Tw = r; , 
6 20 
Da = T’ Dp = T > (D2) 


the transport equations of chemical potentials, i.e., 
left-handed SU(2) doublet tops ji:,2, left-handed SU(2) 
doublet bottoms 4,2, left-handed SU(2) singlet tops 
Ute 2(charge conjugation of tr), Higgs bosons pp.2, are 
given by, 


—3K 4 blih 2 + 3vw Ks bup 2 + 30 po u,2 =0, 
—2K 4 nhh o + Vw Ks nun o + 2u =0, 
QvwK 1 hhh? + Qui, 2 — 3T y (Ht,2 + Hb,2 + 2 Ltc 2 + 2tn,2) 


rhn =0, (D3) 
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30w Kı thia + 30w Ko e(m?2) ut,2 + Buy o — 8D y(Me,2 + Mee,2 + Hh,2) — 6D m (Me,2 + Mee,2) — 80 w (Me,2 — He,2) 
—3P ss[|(1 + OF +) pe,2 + (1 + 91») Mp,2 + (1 — 914) Mee,2] = 3K 7 .0)m? Mi 5 


30w K1 bhp 9 + 3Up 9 — 3I y(Hb,2 + Mee,2 + Maz) — 8D w (Mo,2 — Me,2) — 8P s5[(1 + 941) ht, 


+(14+ 941 »)po2 + (1 — 9K t)ut 2] = 0, 


3vw Ki thre 2 + 30wK2,e(m7)! pte,2 + 3uhe 3 — 3T y(Ht,2 + Ho,2 + 2ee,2 + 2un,2) 


6D m(Ht,2 + Mee,2) 5 


—3D's5[(1 + 9412) Meo + (1 + 9K 1,5) Hoo + (1 — 9K 1,2) oee,2] = 8K 7 8,m2 the 1 , 


-3K44 2 + 30w K5 tU, 2 + 3¥wKo2(m?)'us2 + 3rtut2 = —3vu Kst (m?0,) + 3vwKo20,m3 (m?) 


oe 2Q/,,/ 
—3K 10,27 F,U4 ) 


—3K 4 tbe 2 + 3, Ks tuhe 9 + 30 Ko (Mm?) ute 2 + 3TH ute o = —30y Ks (m? 0Y + 3Uy Ko tbim? (m7) 
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